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Current observational data favor cosmological models which differ from the standard model due 
to the presence of some form of dark energy and, perhaps, by additional contributions to the more 
familiar dark matter. Primordial nucleosynthesis provides a window on the very early evolution 
of the universe and constraints from Big Bang Nucleosynthesis (BBN) can bound the parameters 
of models for dark matter/energy at redshifts of order ten billion. The spectrum of temperature 
fluctuations imprinted on the Cosmic Microwave Background (CMB) radiation opens a completely 
different window on the universe at epochs from redshifts of order ten thousand to nearly the 
present. The CMB anisotropy spectrum provides constraints on new physics which are independent 
of, and complementary to those from BBN. Here we consider three classes of models for the dark 
matter/energy: extra particles which were relativistic during the early evolution of the universe 
("X"); Quintessence models involving a minimally-coupled scalar field ("Q"); models with a non- 
minimally coupled scalar field which modify the strength of gravity during the early evolution of the 
universe ("G"). We constrain the parameters of these models using data from BBN and the CMB 
and identify the allowed regions in their parameter spaces consistent with the more demanding joint 
BBN and CMB constraints. For "X" and "Q" such consistency is relatively easy to find; it is more 
difficult for the "G" models with an inverse power law potential for the scalar field. 



I. INTRODUCTION 

Current and ongoing spaced-based and ground-based 
observational programs have provided increasingly pre- 
cise data enabling us to view the present and recent 
universe with hitherto unprecedented clarity and detail. 
These data have led to a current "standard" model of 
cosmology in which "dark matter" and "dark energy" of 
unknown origins play significant roles. At the same time, 
the earlier, radiation dominated (RD) evolution of the 
universe remains largely hidden from view. These early 
epochs, which may harbor valuable clues to the nature of 
the dark matter/energy are shrouded by the huge optical 
depth of the pre-recombination plasma. As a result, they 
can only be explored indirectly, through the comparisons 
with observations of the predictions of primordial nucle- 
osynthesis ( "big bang nucleosynthesis" : BBN) and of the 
temperature fluctuations in the spectrum of the cosmic 
microwave background radiation (CMB). Studies of BBN 
and the CMB offer a valuable complement to probes of 
the recent structure and evolution of the universe and, in 
conjunction with them, may provide unique constraints 
on competing models for the dark energy. For example, 
although the abundances of the light elements produced 
during BBN depend largely on the universal density of 
baryons, they are also sensitive to the early universe ex- 
pansion rate which, in turn, is determined by the energy 
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density in relativistic particles and, on the strength of 
gravity. Thus, BBN can not only constrain the contribu- 
tion of additional energy density beyond that predicted 
for the standard model of particle physics, it can also 
probe the strength of the gravitational interaction dur- 
ing such early epochs. Similarly, the CMB fluctuation 
spectrum depends not only on the total energy density 
and the magnitude of its relativistic component, but also 
on the expansion rate and the strength of gravity. As a 
result, BBN and the CMB have the potential to distin- 
guish among - or at least constrain - competing models 
for the dark energy some of which leave the strength of 
the gravity invariant while adding to the energy density, 
while others may modify both. 

In this paper we compare and contrast the modifica- 
tions to the standard model BBN and CMB predictions 
in the presence of extra, relativistic energy ("equivalent 
neutrinos"), for those models of "Quintessence" which, 
during RD epochs, contribute a fixed fraction of the rel- 
ativistic energy density, and for non-minimally coupled 
scalar fields (one of whose effects is to alter the strength 
of the gravitational constant G) . After an introduction 
to, and an overview of the non-standard cosmologies ex- 
plored here (§11), the BBN predictions for these three 
general models are compared and contrasted and current 
data is employed to provide constraints on them in §111. 
In §IV an overview is provided of the physical origin of 
the CMB fluctuations in the standard model as a prelude 
to our discussion of the predictions in the non-standard 
cosmologies (§V). The CMB constraints are presented in 
§VI and in §VII they are combined with those from BBN 
to provide joint constraints on the non-standard cosmolo- 
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gics considered here. Our results on the joint BBN-CMB 
constraints arc summarized in §VIII. Unless otherwise 
stated, we use units in which % = c = 8nG = 1. 



II. NON-STANDARD COSMOLOGIES 

Before we begin discussing the implications for BBN 
and the CMB of the "new physics" we are considering it 
is worth spending some time investigating the effects of 
each upon cosmology in general. As a reference we take 
the standard model to be geometrically flat, containing 
three, light neutrinos, baryons along with cold dark mat- 
ter (CDM), and a cosmological constant (A). For nu- 
merical estimates we will often adopt the so-called con- 
cordance values for the density parameters: flM = 0.3, 
Cl\ = 1 — Qm = 0.7, and for the present value of the 
Hubble parameter: Ho = 72 kms _1 Mpc _1 . 



"X" 

In many extensions of the standard models of cosmol- 
ogy and of particle physics there can be "extra" energy 
density contained in new particles or fields, px- Adding 
energy always results in an increase of the expansion rate, 
H, at a given redshift since 




where p is the total energy density. The increased expan- 
sion rate in turn implies that the age of the Universe at 
a given scale factor, 

Jo an 

is younger than in the standard model. If the stress en- 
ergy tensor for X is conserved then the rate of change of 
the energy density for X obeys 

px + 'ZH( Px +P x ) =0 (3) 

where the dot denotes the derivative with respect to cos- 
mic time and Px the pressure. The equation of state, 
Wx = Px/px, is therefore all that is required to deter- 
mine the evolution of px- 

When X behaves like radiation, that is when wx — 
1/3, this additional energy density varies with the scale 
factor, a, as px oc a~ 4 . In the remainder of this paper 
we shall assume this behavior for "A" . Of course, other 
equations of state are possible, leading to different scal- 
ings of the X-density with a but, because this evolution of 
px for w x — 1/3 is exactly the same as that for the pho- 
ton and neutrino energy densities [75] extra relativistic 
energy is a simple extension of the standard model. 



"Q" 

Another example of "new physics" is the positing of 
a Quintessence field to replace the cosmological constant 
as a source of the dark energy in the standard model [1]. 
Quintessence has become an increasingly attractive alter- 
native to A because, by making the dark energy dynamic, 
it helps alleviate the large discrepancy between the mat- 
ter/radiation and vacuum energy densities most appar- 
ent during the early Universe. Virtually all Quintessence 
models are taken to be a minimally-coupled scalar field 
with an energy density pq — <j) 2 /2 + V{4>) and a pressure 
Pq = 4> 2 /2 — V(4>) where V{4>) is the potential energy. 
If the field only interacts with the other constituents of 
the cosmic fluid gravitationally then equation (3) also 
applies to pq but the characteristic of Quintessence is a 
non-trivial equation of state that is not known a priori. 
However by substituting the expression for pq and Pq 
into equation (3) we obtain the Klein-Gordon equation 
for <j> 

4> + 3H<j) + V^ = (4) 

where V$ = dV/d<p. Here, once V{<f>) is specified, there 
are no unknown functions. For many potentials there 
exist solutions to this equation to which the field con- 
verges from a wide range of initial conditions. Solutions 
of this type are generically dubbed "tracker" solutions 
even though this term was originally introduced to dis- 
tinguish a specific class of models [1]. We shall assume 
that the field has reached this tracker solution long before 
BBN. The evolution of the field is controlled by the form 
of the potential and there are many different models in 
the literature from which to choose. This is unfortunate 
because there is then no generic behavior or consequence 
of a dynamic dark energy except that we may state that if 
the Quintessence equation of state wq is ever larger than 
— 1, then the total energy density at a given redshift must 
be larger than in the standard model. If we restrict our- 
selves to tracker potentials then, very roughly, we may 
divide them into those where the Quintessence energy 
density during the radiation dominated epoch is signif- 
icant and those where it is not. Examples of the latter 
type include the popular Inverse Power Law potentials, 
V oc (f)~ a [2, 3, 4, 5]. If a model is of the former type 
then its equation of state must be similar to the equation 
of state of the other components of the cosmic fluid Wf 
and, indeed, it is even possible during some epochs for 
the two to be equal. The Quintessence models we are 
considering are exactly of this type. 

Any potential that fulfills the relationship wq = Wf 
must reduce to the exponential potential [76] [5] V cx 
cxp(— A (j)) though the Inverse Power Law potential fulfills 
this requirement when a — > oo as may other potentials 
in appropriate limits. This popular potential has been 
investigated by Ferreira & Joyce [6], Liddle & Scherrer 
[5] and Copeland, Liddle & Wands [7], amongst others 
who showed that the density parameter of the tracker 
solution for the field [77] is given by = 3 (wf + 1)/A 2 . 



3 



Note the value of is larger during radiation domina- 
tion than during matter domination: CI*q D ^ > O.q 10 ^. 
The increase in energy density results in a swifter expan- 
sion and the age of the Universe at a fixed scale factor, 
again given by equation (2), is smaller than that of the 
standard model by the approximate factor 1 / ^Jl — Qq. 
However, by itself, the exponential potential is incapable 
of leading to an accelerating Universe dominated by dark 
energy if the Quintessence field has reached it's tracker 
solution. In order to account for this observation the 
potential must depart from this simple form as the field 
and the Universe evolve. More precisely, the field must 
become dominated by its potential energy in order for 
wq < —1/3. There are many ways to achieve this re- 
sult by modifying the potential from its pure exponential 
form. One modification, which we will not discuss here, 
was proposed in Dodelson, Kaplinghat & Stewart [9] . In- 
stead, for our analysis we have selected the Albrecht- 
Skordis (AS) model [8] where the potential is of the form 
of a product of a polynomial (quadratic in this case) in 
<j) and the exponential, 



V(4>) = [(</>-0 o ) 2 + A] exp(-A, 



(5) 



The polynomial introduces a local minimum and maxi- 
mum into the potential and the parameter A must satisfy 
A < 1/ A 2 in order for the minimum to be at a real value 
of (f> [10]: for specificity we have chosen AX 2 — 1/100. 
With this choice for A A 2 the two extrema are then very 
close to 0o and 0o + 2/ A. The tracker solution for <ft thus 
evolves according to the pure exponential potential at 
early times when <j> -C 0o. As <p evolves and approaches 
the minimum at </>o the polynomial steepens the poten- 
tial relative to the exponential resulting in an increase 
of (f> and a simultaneous decline in the potential energy. 
This phase of wq > Wf leads to a significant decrease in 
Qq. The field passes through <frg whereupon wq begins to 
decrease as the field climbs out of the minimum. When 
wq = Wf again £Iq reaches its minimum value. If the 
local maximum of the potential is sufficiently high the 
field is unable to pass over the maximum at 4>o + 2/ A and 
resume the exponential scaling behavior, instead it be- 
comes trapped and begins oscillating around (f>o with an 
ever decreasing amplitude. After several oscillations the 
kinetic energy has decreased sufficiently and wq <~ — 1 
at which point the field begins to mimic a cosmological 
constant. 

On a technical note, while it may appear that we have 
two remaining free parameters in V{4>) after specifying 
A A 2 , the solution of equation (4) depends upon H which 
in turn is a function of ^mo & n d Hq. If we fix the ge- 
ometry of the Universe to be flat then we have specified 
the energy density of the field at the present time and so 
introduced a constraint and eliminated another degree of 
freedom. This leaves only one remaining free parameter 
which we choose to be A. Then, assuming the amplitude 
of the oscillations around the minimum are negligible at 
the present time, the imposition of the boundary condi- 



tions leads to the result that X(f>o <x —2 In A. Also, from 
our numerical calculations, we find that the minimum 
value of Qq is proportional to A 2 which implies that the 
rcdshift at which the field begins to resemble a cosmolog- 
ical constant, z\, increases with A. We find that roughly 
z A ss 4.7 A 2 / 3 or equivalently z A w 7.5 ( Q ( q D) )- 1 / 3 . 



A third example of "new physics" is the more radical 
proposal of a non-minimal coupling between the Ricci 
scalar, R, and a scalar field which we will further pro- 
mote to the role of a Quintessence dark energy. This 
extension therefore requires the postulation of the form 
of the coupling in addition to the potential V{4>). Again 
there are many different models in the literature from 
which to choose each involving a different form for the 
coupling and the potential. A very interesting general 
feature of this class of models is that the strength of the 
effective gravitational constant during the early evolution 
of the universe may differ from its current value, which 
is fixed by terrestrial and solar system experiments. For 
this reason, we label these models by "G". While it is 
possible to derive a general formalism for the cosmology 
and the evolution of perturbations [11, 12, 13] there is 
no universal behavior and therefore it becomes necessary 
to restrict ourselves to a specific example. The model we 
have adopted here is the minimal extension of the Non- 
Minimally Coupled (NMC) model investigated by Chen, 
Scherrer & Steigman [14] and Baccigalupi, Matarrese & 
Perrotta [15] amongst others. In this model the action 
takes the form 



d 4 x- 



F(4>)R 



V(<f>)+C f 



(6) 



with F{4>) = 1 + £(</> 2 — </>o)i 0o is the value of the field 
today, and £ is the coupling constant. From the action 
we can define the cosmological gravitational parameter 
to be 1 /F and so the evolution of the field, and therefore 
F, will lead to an evolution of the strength of gravity. 
The potential is taken to be the previously mentioned 
Inverse Power Law V — Vo 4>~ a which is known to be a 
viable Quintessence model in the minimally coupled limit 
[5]. There are two distinct approaches to modifications 
that arise in these scenarios: redefine the energy density 
and pressure of the field leaving the cosmological equa- 
tions unaltered or, adopt the minimally coupled defini- 
tions and work with the modified cosmological equations. 
Both approaches are, of course, equivalent and here we 
adopt the latter. As a result, the Friedmann equation is 
modified, becoming 



H 2 + 



HF 



P 
3F 



and, after introducing the function E, defined by 

3F 2 

E= 1 + 



2F 



(7) 



(8) 
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the field evolves according to the equation 

Fa R 



+ 3i?( 

Fa 



Vcb = 



2FE 



( P -3P + 3FaVa -3F 00 2 ) . (9) 



The presence of the right hand side of equation (9) has 
been dubbed the R-boost [15]. Examining the right hand 
side of equation (9) we discover that the last two terms in 
parentheses on the right-hand side of the equation com- 
bine to give a contribution of 



F A V A 
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-£(l + a)p G -£(l-a)P G (10) 



for this potential and coupling. Again we must enforce a 
self-consistent cosmology because the Hubble parameter 
is still a function of I^mo and H so we must adjust the 
normalization constant in the potential, Vo, to ensure the 
boundary conditions are matched. 

The Inverse Power Law potential in the minimally cou- 
pled limit (£ — > 0) is a well known and frequently studied 
Quintessence model [5]. The energy density of the field 
is negligibly small during the early Universe so that H 
is dominated by the radiation and matter densities until 
close to the present time. The tracker solution for (f> is 
4> oc a 3 ( M, /+ 1 )/( Q! + 2 ) and for a = (9(1) we can immediately 
see that, for the tracker solution, (p <C tfio during much of 
the evolution of the Universe. The equation of state for 
the field is 



w G 



oeWf 



(11) 



which is always smaller than Wf and so the Quintessence 
energy density grows relative to the matter + radi- 
ation fluid [78]. This solution is also stable in the 
sense that perturbations from the tracker behavior are 
damped [5]. The energy density of the field scales as 
pa oc a - 3a ("'/+ 1 )/( a + 2 ). As we approach the present 
era the equation of state for the field begins to deviate 
from the analytic formula in equation (11) and begins 
to approach wq = —1- The emergence of pa, together 
with the descending wq, launches a phase of cosmologi- 
cal acceleration. As a increases the redshift at which the 
departure from equation (11) occurs also increases but 
the redshift at which the Universe begins to accelerate 
falls because while the Quintessence equation of state 
is no longer given by equation (11) its present value is 
still correlated with a with smaller a leading to values 
of wq that are closer to —1 [19]. For the concordance 
model parameters f^Mo = 0.3, flco = 1 — ^mo = 0.7, 
Hq = 72 kms _1 Mpc _1 , the equation of state at the 
present time is larger than —1/3 for a > 8 while the 
requirement that the Universe be accelerating requires 
a<4[19]. 

As with our minimally coupled case "Q" we shall as- 
sume that the field has reached its tracker solution long 
before BBN and so our initial condition for (f> is this limit 
during the radiation dominated epoch. This may seem 



a minor point but the evolution of the field in its non- 
tracker state will be very different from its behavior in 
the tracker solution. We posit that the tracker solutions 
for non-zero £ are the same as in the minimally coupled 
limit and so we must show that the R-boost is negli- 
gibly small at early times (i.e., it diverges with a at a 
slower rate than Va as a — > 0). The coupling F is al- 
most constant and Fa oc cj) is very small so the scaling of 
the R-boost during this phase is controlled by whether 
Pm oc a~ 3 or pa oc a~ 4a /("+ 2 ) is the more divergent: for 
a < 6 it is pm otherwise it is pa- So we see that the two 
terms in equation (9) scale as Va oc a~ 4 ( Q+1 )/( Q+2 ) while 
Fa Rot a -( 2 +3«)/(«+ 2 ) for a < 6, oc a (4-4a)/(a+2) 

for a > 6. Thus the derivative of the potential diverges 
more rapidly than the R-boost whatever the value of a 
and the tracker solution during radiation domination is 
identical to the solution in the minimally coupled limit. 
It could be argued that pu — ► as a — > because there 
are no non-relativistic particles at such high tempera- 
tures in which case R oc a~ 4a / ( a + 2 ) independent of a but 
this makes no difference to our conclusion. This tracker 
solution is also stable for any value of £ as shown by Bac- 
cigalupi, Matarrese and Perrotta [15] who relied on the 
fact that F is virtually constant in order to generalize the 
result of Uzan [18]. During matter domination the situ- 
ation is slightly different: <j) oc a 3 ^ a+2 ^ so both Va and 
FaR scale as a~ 3 ( Q+1 '/( Q+2 ) and, therefore, the tracker 
solution is again the minimally coupled tracker but the 
normalization changes to account for the presence of the 
R-boost. 

Since Fa R scales more slowly than the derivative of 
the potential during the radiation dominated epoch it 
might be expected that eventually an R-boost phase will 
occur: we can estimate the scale factor at which the two 
are equal to be a <~ a eq /^ where a eq is the scale factor at 
radiation-matter equality. An R-boost phase will only 
occur if £ > 1 and, once initiated, will continue into 
matter domination. If £ is small then no R-boost phase 
will occur during radiation domination and neither will it 
commence during the matter dominated epoch. R-boost 
phases can occur when the initial value of <j) differs from 
the tracker value. In this case, if we again assume F is 
almost constant and again approximate i? as R = pm / F, 
the solution for 4> during the radiation dominated epoch 



is 



h 2/3 



for £ < 0, (12) 



for £ > (13) 



where (3 2 = 3 |£| fijw* and 0* is the value of the field at a*, 
Ji and I\ are the Bessel and modified Bcssel functions of 
the first with index 1. In either case, when the argument 
of the Bessel function is small 6 evolves as 



2 W* 



(14) 
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the minus (plus) sign is for £ < (£ > 0). Note that 
when £ < the field moves backwards. The very slow 
change of <j> justifies the assumption that F is almost 
constant but, in contrast with the tracker solution, here 
the behavior arises because <j> is essentially fixed rather 
4> <C <j)Q- During this R-boost phase 4> oc 1/a and the 
potential energy is much smaller than the kinetic energy 
so that wq — > 1 but pa oc 1/a 2 , very different from the 
minimally coupled limit. Eventually any R-boost phase 
will terminate and the field will follow the tracker solution 
but, if 4>+ is much larger than the tracker value, there may 
be a very long delay before this occurs. The approximate 
scale factor at which the R-boost phase ends is given by 
a = a* (^/^(a,))(«+ 2 )/ 4 . 

With the coupling F almost constant during the early 
Universe and the Quintessence energy density entirely 
negligible equation (7) essentially reduces to the Fried- 
mann equation of our standard model except for an ef- 
fective gravitational strength G' = G/F. Therefore the 
only change is to the age of the Universe at a given scale 
factor, once again given by equation (2), which is simply 
rescaled by the factor yfF. If F were not constant, and 
could be adequately described by a power law function of 
the scale factor, then we would attain the circumstances 
investigated by Carroll and Kaplinghat [17]. 

The evolution of the field becomes more complicated as 
the present epoch is approached: the R-boost is no longer 
negligible, F starts to change noticeably and the field's 
energy density becomes important. Again the equation 
of state for the field begins to descend towards wq = — 1 
but the exact evolution is now also a function of £. This 
can be understood from the evolution of pc which is 



p G + 3H(p G + P G ) = ™. 



(15) 



The power source F R/2 is proportional to £ so if £ > 
then the field gains energy relative to a field in the 
minimally coupled limit if R > 0. The power source 
will become a drain if or when R switches sign as the 
Universe begins to accelerate. The increase in energy 
density results in an increase in the equation of state 
relative to the minimal case and thus a relative increase 
in (j> and <j>. The increased energy density will terminate 
matter domination at a higher redshift. 

Finally we impose the constraints on the model param- 
eters from the timing experiment using the Viking probe 
and from limits to the evolution of the strength of gravity 
[16]. The first constraint, the more severe of the two, is 



£</>o < 0.022 
while the second limits 



2£« 



»o (pa 



< 10 



-li 



yr 



(16) 



(17) 



where <po is time derivative of the field at the present 
epoch. As shown in Chen, Scherrer and Steigman [14], (£>o 
is very weakly dependent upon £ so that it is essentially 



determined by a and Vo- Therefore, with fixed values 
of a and Vo, the limit in equation (16) is essentially a 
limit on £ or, if Fbbn is the (almost constant) value of 
F during BBN and if <f> during this epoch is much smaller 
than <po, then 



1 — Fbbn 



(18) 



and the limit on £ therefore becomes a limit on Fbbn or, 
equivalently, on G'/G. From our numerical calculations, 
again with the concordance model parameters, we find 
that G'/G is restricted to lie between 



0.976 < G'/G < 1.025 a = 1, 

0.964 < G'/G < 1.040 a = 2, 

0.942 < G'/G < 1.067 a = 4. 



(19) 



III. 



BIG BANG NUCLEOSYNTHESIS 



To better appreciate the similarities and differ- 
ences among the three candidates for "new" physics 
under consideration here (equivalent neutrinos: "A"; 
Quintessence: "Q"; non-minimal coupling: "G") it will 
be helpful to briefly review "standard" BBN (SBBN). To 
this end the discussion may begin when the universe is 
a few tenths of a second old and the temperature is a 
few MeV. The energy density receives its dominant con- 
tributions from the relativistic particles present; prior to 
annihilation these are: CBR photons, e ± pairs, and 
three flavors of neutrinos, 
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PR = P~f+ Pe+ 5pv = -rP~f- 



(20) 



At this time (T ~ few MeV) the neutrinos are beginning 



to decouple from the photon 



plasma and the neu- 



tron to proton ratio, crucial for the primordial abundance 
of 4 Hc, is decreasing. As the temperature drops below ~ 
2 MeV, the two-body collisions interconverting neutrons 
and protons become too slow to maintain equilibrium and 
the neutron-to-proton ratio begins to deviate from (ex- 
ceeds) its equilibrium value ([n/p) eq = exp(— Am/T)). 
Prior to e ± annihilation, at T m 0.8 MeV when the uni- 
verse is ~ 1 second old, the two-body reactions regulating 
the n/p ratio become too slow compared to the univer- 
sal expansion rate and this ratio "freezes in", although 
it actually continues to decrease due to the emerging im- 
portance of ordinary beta decay (r„ = 885.7 sec). Since 
there are several billion CBR photons for every nucleon 
(baryon), no complex nuclei exist at these early times. 

BBN begins in earnest after annihilation, at T rts 
0.08 MeV (t ~ 3 minutes), when the number density 
of CMB photons with enough energy to photodissociate 
deuterium (those in the tail of the black body distribu- 
tion) is comparable to the baryon density. By this time 
the n/p ratio has further decreased due to beta decay, 
limiting (mainly) the amount of helium-4 which can be 
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synthesized. As a result, the predictions of primordial 
nucleosynthesis depend sensitively on the early expan- 
sion rate. In SBBN it is assumed that the neutrinos are 
fully decoupled prior to e ± annihilation and don't share 
in the energy transferred from the annihilating e ± pairs 
to the CMB photons. Thus, in the post-e 1 * 1 annihilation 
universe the photons are hotter than the neutrinos and, 



PR = P~< + 3p„ = 1.6813p 7 



(21) 



During these RD epochs the age and the energy density 
are related by |pR,i 2 = 1 (recall that we have chosen 
units in which 8nG = 1), so that the age of the universe 
is known (as a function of the CMB temperature) once 
the particle content (pr) is specified and the strength of 
the gravitational interaction (G) is fixed. In the standard 
model, 



Pre — e annihilation 



t T 2 = 0.738 McV 2 s, (22) 



Post - c ± annihilation : t T 2 = 1.32 MeV 2 s. (23) 

The BBN-predicted abundances of deuterium, helium-3 
and lithium are determined by the competition between 
various two-body production/destruction rates and the 
universal expansion rate, while the helium-4 abundance 
depends most directly on the neutron abundance at the 
time BBN begins. As a result, the D, 3 He, and Li abun- 
dances are sensitive to the post-e ± annihilation expan- 
sion rate, while that of 4 He depends on both the pre- 
and post-e ± annihilation expansion rates; the former de- 
termines the "freeze-in" and the latter the importance of 
beta decay. Of course, the BBN abundances do depend 
on the baryon density (7710 = 10 10 n B /n 7 = 274£l B /i 2 ), so 
that the abundances of at least two different relic nuclei 
are needed to break the degeneracy between the baryon 
density and a possible non-standard expansion rate re- 
sulting from new physics or cosmology. 



A. Non-Standard BBN 

Our simplest alternative to the standard cosmology is 
the scenario of extra relativistic energy denoted by "A" . 
When A is decoupled in the sense that it doesn't share 
in the energy released in e ± annihilation, it is convenient 
to account for this extra contribution to the standard- 
model energy density by normalizing it to that of an ex- 
tra, "equivalent" neutrino flavor [20], 



7 

px = AN vPv = -AA„p 7 . 



(24) 



particle-antiparticle pairs annihilated (or unstable parti- 
cles decayed). In this case, the contribution to AN U from 
each such particle will be < 1 (< 4/7). Since we are in- 
terested in different models of non-standard physics re- 
sulting in modifications to the standard model energy 
density and expansion rate, henceforth this case will be 
identified by a superscript, A; AAf„ = AN* . 

In the presence of this extra component, the pre- 
e ± annihilation energy density in equation (20) is modi- 
fied to, 



(PR)pr 



43 



pre 



7AN* \ 
43 J 



p 7 . 



(25) 



The extra energy density speeds up the expansion of 
the universe so that the right hand side of the time- 
temperature relation in equation (22) is smaller by the 
square root of the factor in parentheses in equation (25). 
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-'pre 



(t/f) 



1 



7 AN, 



x\ 1/2 



1 ■ 1 43 
= (1 + 0.1628AA*) 1 / 2 . (26) 



In the post-e^ 1 annihilation universe the extra energy 
density contributed by the As is diluted by the heating 
of the photons, so that, 



( PR )* = 1.6813 (1 + 0.1351AAT* )p 7 



and 



^post 



(t/t') post = (1 + 0.1351AA*) 1 



/2 



(27) 



(28) 



These relations (eqs. 25-28) may now be generalized to 
the two other cases under consideration. 

For our minimally coupled Quintessence model "Q" 
the energy density of the field during radiation domina- 
tion (relevant to BBN), is fl ( g D) = 4/A 2 . This extra 
energy density may be written in terms of an equivalent 
A A,?, 



AJV? = « 
7 



n 



(RD) 



i-n { Q RD \ 



(29) 



Since, for this Quintessence model, the scalar field con- 
tributes the same fraction of the total (radiation) energy 
density pre- and post-e ± annihilation, 



/pre 



(PR)po 



43 



1 



7ANQ\ 

- j Pn 



43 



(30) 



and 



For each such "neutrino-like" particle (i.e., a two- 
component fermion), if Tx = T v , then AA„ = 1; if A is 
a scalar, AA„ = 4/7. However, it may well be that A 
has decoupled earlier in the evolution of the universe and 
has failed to profit from the heating when various other 



S£ e = S® st ^SQ = (1 + 0.1628AA?) 1 



/2 



(31) 



Thus for this class of Quintessence models the speed- 
up factor (5"2) in the universal expansion rate prior 
to e ± annihilation is the same as the speed-up factor 
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after e annihilation. In comparing with the equiv- 
alent neutrino case, we see that for AA„ = AN? — 
AN?, the post-e 1 * 1 annihilation universe expands faster 
for "Q" than for "A"; alternately, for the same post- 
e ± annihilation speed-up, AN® « 0.83AA* . 

In our non-minimally coupled Quintessence model, 
"G" , the Quintessence energy density during BBN is en- 
tirely negligible hence the total energy density at a given 
redshift is unaltered, 

(PR.) G - PR (32) 

where /3r is given by eqs. 25 and 27 for the pre- and post- 
e ± annihilation universe respectively. However, since t oc 
(G P n)- 1/2 , 

S° e = S? ost = (G'/G)^ ee (l + ^) V2 . (33) 

For this case we have defined an equivalent number of 
extra neutrinos by 

. Arr 43 / AG\ 

AN? ee — i^—j . (34) 

As with Quintessence, the non-minimally coupled 
scalar fields result in faster expansion of the post- 
ed 1 annihilation universe for the same increase in the 
pre-e^ annihilation expansion rate. Thus, for the same 
AN V the effects on BBN of "Q" and "G" are identical, 
but they do have the potential to be different from those 
due to the usual example of extra energy density in the 
form of equivalent neutrinos ("X"). However there is one 
important difference between "Q" and "G" , namely that 
AN? is unconstrained (at the moment) while the limits 
expressed in (16) and (17) translate to a restricted range 
for AN?. The range for AN? is a function of a, the ex- 
ponent of the Inverse Power Law potential and, as shown 
in Chen, Scherrer and Steigman [14], the limits increase 
with the exponent. 

On the basis of the above discussion the results of our 
detailed BBN calculations may be understood. Since the 
primordial abundances of D, 3 He, and Li freeze-in late, 
well after annihilation has occurred, they mainly pro- 
vide a probe of S pos t- In contrast, the He mass fraction, 
Yp, retains sensitivity to both S pre and S pos t- Further- 
more, while the abundances of D, 3 He, and Li are most 
sensitive to the baryon density, the 4 He mass fraction 
provides the best probe of the expansion rate. This is 
illustrated in Figure 1 where, in the AN V - r/io plane 
are shown isoabundance contours for D/H and Yp (the 
isoabundance curves for 3 He/H and for Li/H, omitted for 
clarity, are similar in behavior to that of D/H). While 
the solid curves are for the "usual" extra energy density 
in relativistic particles, the dashed curves show, for the 
same equivalent AN V (see eqs. 29 & 34), the modifica- 
tions to SBBN for both minimally and non-minimally 
coupled scalar fields. The trends illustrated in Figure 1 
are easy to understand in the context of our discussion 
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FIG. 1: Isoabundance curves for D and 4 He in the AN V - 
7710 plane. The solid curves are for AN„ = AN*, while 
the dashed curves correspond to the minimally and non- 
minimally coupled scalar field cases where AN V = AN? = 
AN?. For 4 He the curves (nearly horizontal) are for, from 
top to bottom: Y = 0.25, 0.24, 0.23. For D the curves (nearly 
vertical) are for, from left to right: 10 5 (D/H) = 4.0, 3.0, 2.0. 



above of SBBN. The higher the baryon density (r?io), 
the faster primordial D is destroyed, so D/H is anticor- 
related with r]i . But, the faster the universe expands 
(AN V > 0), the less time is available for D-destruction, 
so D/H is positively correlated with AN U . Since for the 
same AN V , the post-e 1 * 1 annihilation universe expands 
more rapidly for «<j>" and "G" than for "A" , the D- 
isoabundance curves differ slightly as shown in Figure 1. 
In contrast to D (and to 3 He and Li), the 4 He mass frac- 
tion is relatively insensitive to the baryon density, but is 
very sensitive to both the pre- and post-e 1 *" annihilation 
expansion rates (which control the neutron-to-proton ra- 
tio) . The faster the universe expands, the more neutrons 
are available for 4 He. Again, the effect on Yp of the same 
AN U is slightly different for "A", than for "Q" and "G", 
as seen in Figure 1. 



B. Primordial Abundances 

It is clear from Figure 1 that any BBN constraints 
on new physics will be data-driven. While D (and/or 
3 He and/or Li) largely constrain the baryon density and 
4 He plays the same role for AN V , there is an interplay 
between 7710 and AN V which is quite sensitive to the 
adopted abundances. For example, a lower primordial 
D/H increases the BBN- inferred value of 7710, leading 
to a higher predicted primordial 4 He mass fraction. If 
the primordial 4 He mass fraction derived from the data 
is "low", then a low upper bound on AN„ will be in- 
ferred. It is, therefore, crucial that we make every effort 
to avoid biasing our conclusions by underestimating the 
uncertainties at present in the primordial abundances de- 
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rived from the observational data. To this end, first of 
all we concentrate on deuterium as the baryometer of 
choice since its observed abundance should have only de- 
creased since BBN [21] and the deuterium observed in the 
high rcdshift, low mctallicity QSO absorption line sys- 
tems (QSOALS) should be very nearly primordial. The 
post-BBN evolution of 3 He and of Li are likely more com- 
plicated, involving competition between production, de- 
struction, and survival. 

Even so, inferring the primordial D abundance from 
the QSOALS has not been without its difficulties, with 
some abundance claims withdrawn or revised. At present 
there are 4-5 QSOALS with reasonably firm deuterium 
detections [22, 23, 24, 25, 26]. However, when D/H is 
plotted as a function of metallicity or redshift, there is 
significant dispersion and the data fail to reveal the an- 
ticipated "deuterium plateau" [27]. Furthermore, subse- 
quent observations of the D'Odorico et al. [26] QSOALS 
by Levshakov et al. [28] revealed a more complex veloc- 
ity structure and led to a revised - and uncertain - deu- 
terium abundance. This sensitivity to the often poorly 
constrained velocity structure in the absorbers is also ex- 
posed by the analyses of published QSOALS data by Lev- 
shakov and collaborators [29, 30, 31], which lead to con- 
sistent but somewhat higher deuterium abundances than 
those inferred from "standard" data reduction analyses. 
Given this current state of affairs we believe that while 
the O'Meara et al. [24] estimate for the primordial abun- 
dance is likely accurate (D/H = 3.0 x 10~ 5 ), their error 
estimate (±0.4 x 10~ 5 ) may be too small. At the risk of 
erring on the side of caution, here we adopt a range which 
encompasses the uncertainties referred to above: D/H = 
3.0+ ' 5 x 10~ 5 . Although permitting a larger than usual 
range in baryon density, this choice has little direct effect 
on the probe of new physics (constraints on AN U ) which 
is the focus of this study. 

A similarly clouded situation exists for the primordial 
abundance of 4 He. At present there are two estimates for 
the primordial abundance of 4 He based on large, nearly 
independent data sets and analyses of low-metallicity, 
extragalactic Hn regions. The "IT" [32, 33] estimate 
of Yp(IT) = 0.244 ± 0.002 and the "OS" determina- 
tion [34, 35, 36] of Yp(OS) = 0.234 ± 0.003 which is 
nearly 3a lower. Recent high quality observations of a 
relatively metal-rich Hn region in the Small Magellanic 
Cloud (SMC) by Peimbert, Peimbert, and Ruiz (PPR) 
[37] reveal an abundance Y SM c = 0.2405 ±0.0018. When 
this abundance is extrapolated to zero metallicity, PPR 
find Yp(PPR) = 0.2345 ± 0.0026, lending some support 
to the lower OS value. These comparisons among dif- 
ferent observations suggest that unaccounted systematic 
errors may dominate the statistical uncertainties. In- 
deed, Gruenwald, Steigman, and Viegas [38] argue that 
if unseen neutral hydrogen in the ionized helium re- 
gion of the observed Hn regions is accounted for, the 
IT estimate of the primordial abundance should be re- 
duced to Yp(GSV) = 0.238 ± 0.003 (see also [39, 40]). 
Here, we adopt this latter estimate for the central value 



but, as we did with deuterium, the uncertainty is in- 
creased in an attempt to account for likely systematic 
errors: Y P = 0.238 ± 0.005, leading to a 95% CL range, 
0.228 < Yp < 0.248; this is in agreement with the es- 
timate adopted by Olive, Steigman, and Walker [41] in 
their review of SBBN. 



C. Constraints From BBN 
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FIG. 2: Likelihood contours (68%, 95% and 99% respectively) 
in the AN V - 7710 plane. The solid curves are for AJV„ = 
AN X , while the dashed curves correspond to the minimally 
and non-minimally coupled scalar field cases where AN V = 
AN^f = AN®. The best fit points are indicated by the cross 
for "X", and the square for "Q"/"G". 



Using the abundance ranges adopted above, we have 
calculated the AA„ - likelihood contours which are 
shown in Figures 2 & 3. A simple, semi-quantitative 
analysis can serve to shed light on the detailed results 
shown in Figures 2 & 3. As revealed in Figure 1, the 
baryon density is primarily fixed by deuterium. For D/H 
= 3.0 x 10" 5 , 7710 « 5.6 (Q B w 0.020). In contrast, AN„ is 
most sensitive to 4 Hc; AY« 0.013AA„. For 7710 ~ 5.6 
and AN U = 0, Yp = 0.247. Comparing this with the cen- 
tral value (0.238) or the 2a upper bound (0.248) inferred 
from the data, it can be expected that AA„ w —0.7 and 
AN U < +0.1 respectively. In fact, the detailed calcu- 
lations reveal that the "best" value for "X" is AA„ = 
-0.65, while for "G/Q", AN U = -0.58; the 2a upper 
bounds are AN U < +0.04 in all three cases, correspond- 
ing to the bound ^q ?D - ) < 0.007 for the Quintessence 
model (A > 25). If in place of our adopted central 
value and range for Yp, those from IT had been cho- 
sen, the contours in Figures 2 & 3 would shift upwards by 
AA„ w 0.5 and they would be narrower in the vertical di- 
rection. However, since at 95% CL the two estimates for 
Yp agree, the upper bounds to AA„ will be closely equal. 
While we are most concerned with the BBN constraints 



9 



on AN V , we note in passing that the best estimate for 
the baryon density is r/ w w 5.2 (Q B h 2 w 0.019). 

1 | , , 1 , 1 , 1 , 1 , 1 1 



epoch. We turn next to the effects on the formation of 
the CMB anisotropy spectrum and the constraints which 
may be imposed on these models by the current observa- 
tional data. 




FIG. 3: Likelihood contours as in Figure 2, but with the 
restriction that AN„ > 0. The best fit points are indicated 
by the cross for "X" , and the square for "Q" / "G" . 



We must point out that there is a logical inconsistency 
in the above analysis. Firstly, for two of the cases we (and 
others before us) have been considering, "A" and "Q", 
the energy density and expansion rate during BBN are 
only increased compared to those for SBBN [79] . Since it 
is well established that there are three, very light (hence, 
relativistic at BBN) neutrinos, N v > 3 and A7V„ > 0. 
Only our "G" case can naturally accommodate S < 1 
and, hence, AN V < (see, e.g., Chen, Scherrer, and 
Stcigman [14]). In Figure 3 are displayed the correspond- 
ing likelihood contours when the restriction AN V > 
is imposed. The "best fit" now occurs for AN V = 
and r/io = 5.0 {VL^h 2 = 0.019). In this case the up- 
per bounds on AN„ are increased relative to those when 
AN V is left free to be negative as well as positive; at 
95% CL, AA„ < 0.43 for "A" and AN„ < 0.39 for 
"Q". For the Quintessence model this corresponds to 
the bound CI { q D) < 0.060 (A > 8.2). From a similar 
analysis (AN V > 0), but using the narrower range of the 
IT 4 He abundance as well as the O'Meara et al. [24] nar- 
row deuterium range, Bean, Hansen, and Melchiorri [44] 
find a 95% upper bound of Q Q < 0.045 (A > 9.4). 

Secondly, as mentioned before in §11, the range of AN^ 
is constrained by the experimental limits on the present- 
day deviation from GR as indicated by the allowed ranges 
of G listed in (19). The permitted ranges of AA,^ are, 
therefore, a function of a. Without (yet) imposing a 
restriction on this parameter, the range in AN V is other- 
wise unconstrained. 

As the results here demonstrate, BBN can impose sig- 
nificant constraints on these non-standard models, no- 
ticeably restricting the choice of the additional param- 
eters which accompany the new physics. The changes 
in the dynamics and evolution of the universe imposed 
by these models continues after BBN, up to the present 



IV. STANDARD CMB 

In order to understand the effects of new physics upon 
the formation of the anisotropies in the CMB we briefly 
review the formation of temperature fluctuations in the 
standard model. More extensive reviews may be found 
in Ma and Bertschinger [45], Hu and Sugiyama [46] and 
elsewhere so we will only outline the general procedure 
and important results. Hereafter overdots denote deriva- 
tives with respect to conformal time r where a(r) dr = dt 
and overbars indicate unperturbed quantities. We shall 
write the conformal Hubble parameter as H so that 
aTL = da/dr and therefore 

r{a)= I do. (35) 
Jo a H 

The CMB anisotropies observed today are the red- 
shifted temperature fluctuations that occurred in the 
baryon-photon fluid prior to recombination. The stress 
energy for a perfect fluid with energy density p and pres- 
sure P is simply = P g£ + (p + P) u v and the 
perturbations in the fluid are introduced asp = p + 5p = 
p(l + 5), P = P + 5P and = + but also an 
anisotropic stress £ M!y = PIP" is included in order to 
account for shear, viscosity and other such processes [45]. 
For multiple components the stress energy tensor is the 
sum of the stress energy tensors of each component i, in 
our case vacuum (A), CDM (C), baryons (B), photons 
(7) and neutrinos (y), so that p — X^Pi, etc. For the 
relativistic components, such as photons and neutrinos, 
the treatment as a fluid is unrealistic and, instead, the 
phase space description must be used. The temperature 
fluctuation, or brightness function, 6^ = STi/Ti, is intro- 
duced for each relativistic species; the &i are functions of 
time r, position x and the direction cosines n of the par- 
ticle momenta. There are no fluctuations in the vacuum 
energy. 

The manner in which scalar perturbations in the met- 
ric are introduced defines the gauge and many differ- 
ent choices have appeared in the literature. Two of the 
most popular gauges, the synchronous and the Newto- 
nian, were compared in Ma and Bertschinger [45] and 
others are discussed by many authors including Hwang 
[47]. The synchronous gauge is chosen so that the pecu- 
liar velocity of the CDM is zero. 

The perturbations in all quantities may be decomposed 
into their Fourier modes (with wavenumber A:) and, in 
addition, the angular dependence in the brightness func- 
tions may expanded as an infinite sum of Legendre poly- 
nomials with coefficients 0£. In terms of these moments, 
the thermodynamic equivalents of the energy density, 
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pressure, and anisotropic stress are 



5pi = 4 <3 l0 Vi = On 



12 

n 2 = — e, 2 

5 



(36) 



for each relativistic component i. 

All that remains is to specify the initial conditions for 
the perturbations: this is accomplished by examining the 
behavior on superhorizon scales deep within the radiation 
dominated epoch, retaining only the growing modes, and 
relating the perturbations with some 'principle' such as 
adiabatic or isocurvature. For adiabatic initial conditions 
in the synchronous gauge we have the well known rela- 
tions [45]: 



Sc = 



5 B 
4 



Pb vb 

n„ 



3 A 3 

4 7 4 
15 + 4/,, 



C(fcr) 2 



3 P ^ 7 23 + 4/, 
-16<5c 



(37) 



3(15 + 4/,) 



where /„ is the fraction of the relativistic energy density 
in neutrinos and C is a dimensionless constant. 

The evolution of each mode is a function of its wave- 
length and of the gauge. All modes begin at superhori- 
zon scale, fcr -C 1, and the CDM perturbations evolve 
as Sc oc k 2 r 2 during both matter and radiation domina- 
tion. As time progresses the horizon becomes bigger than 
the wavelengths of different modes; modes with larger 
wavenumbers (smaller wavelengths) enter the horizon be- 
fore the smaller wavenumber (larger wavelength) modes. 
The evolution of Sc for modes that enter the horizon dur- 
ing radiation domination is stunted; they grow only as 
lnr. This stagnation continues until matter domination 
whereupon the modes begin to grow again as r 2 . Prior 
to recombination the behavior of the baryon-photon fluid 
changes radically as a mode enters the sound horizon s 7 #, 
which is given by 



S~fB{T) = dTCj B = -^r; 




Pi 



p~ ( + P 1 + PB 



(38) 



where c 7 s is the baryon-photon fluid sound speed. Note 
the presence of the baryon density in the denominator 
of equation (38). The perturbations present in the fluid 
drive oscillations in the baryon density, the peculiar ve- 
locity, and the photon brightness function for moments 
I € {0, 1}. There are no fluctuations for I > 2 because 
of the tight coupling between the photons and baryons 
and because the isotropising effect of Thompson scat- 
tering suppresses these moments. The oscillations are 
forced by the CDM perturbations via the metric poten- 
tials and their amplitude increases with the ratio of the 
baryon and photon densities. Once initiated, the oscil- 
lations continue until recombination at which point the 
photon-baryon/electron interaction ceases and the pho- 
tons free-stream. The phase of each mode at recombina- 
tion is recorded in the anisotropy angular power spectrum 



Cg. Recombination, which occurs at t+, a*, introduces 
an important modification to the Cg known as Silk damp- 
ing [48] that is due to the increase in photon mean free 
path as the number of free electrons crashes. The sub- 
sequent increase in photon diffusion smoothes temper- 
ature fluctuations on scales smaller than this diffusion 
length. Modification of the 9^ does not end at recombi- 
nation: the evolution of the Universe after this time will 
imprint itself upon brightness function moments through 
the evolution of the potentials. This change, known as 
the Integrated Sachs- Wolfe (ISW) effect, is initiated by 
the fading importance of the radiation density (the early 
ISW effect) and the growing importance of the vacuum 
energy towards the present epoch (the late ISW effect). 
Both cause an increase in the photon temperature vari- 
ance on the scales at which these processes occur. 

What emerges from the detailed calculations is the so- 
lution for the evolution of the perturbations, for a given 
cosmology. The perturbations for each wavenumber are 
scaled according to an initial power spectrum and nor- 
malized, usually to COBE. The photon temperature vari- 
ance vector Cg, or Cg = £(£+ l)Cg/2w, is then con- 
structed. The spectrum of the Cg consists of a series of 
peaks at specific t and so the gross features of a model 
can be characterized by the positions of the peaks, the 
height of the first peak relative to the COBE normal- 
ization point and then the relative heights of all other 
peaks to the first [49, 50]. Regardless of the curvature 
of the Universe, the position of the peaks is proportional 
to the ratio of the comoving angular diameter distance 
and the size of the sound horizon at last scattering [49] . 
The height of the first peak is a measure of change in the 
Universe since last scattering, while the relative heights 
of the other peaks and the ratios of their separations is 
set by the astrophysics prior to and during recombina- 
tion. In general there is a complex interplay between the 
different elements contributing to the formation of the 
temperature fluctuation spectrum, but whenever possi- 
ble we shall try to couch the changes to the CMB in this 
language. 



NON-STANDARD CMB 



•'X r 



The inclusion of extra relativistic energy into the CMB 
calculation is relatively straightforward. As for the pho- 
tons and the three known neutrino flavors, a brightness 
function is introduced that is Fourier transformed and 
the angular dependence expanded in a basis of Legcndrc 
polynomials. If we do not allow the extra energy density 
to interact with the other components (except through 
gravity) then the set of equations governing the evolution 
of the brightness function coefficients is exactly the same 
as for the neutrinos and there is no suppression of the 
higher moments of the brightness function. The initial 
conditions for the extra energy density are taken to be 
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exactly those of the neutrinos so that N v — > N v + AN* 
and /„ — > f v + fx ■ The effects of extra relativistic energy 
density "A" upon the CMB anisotropy power spectrum 
are then equivalent to an amplification of the neutrino 
sector. The extra energy will increase the expansion 
rate H, with the largest changes, relative to the stan- 
dard model, occurring during the early, radiation dom- 
inated, epoch. The swifter expansion modifies both the 
time-temperature relationship and leads to an earlier de- 
coupling of the photons since the scattering rate now be- 
comes smaller than H at an earlier epoch. Therefore 
both age of the Universe, r+, and the scale factor, a*, at 
recombination are smaller and from equation (38) we see 
that this will lead to a smaller sound horizon, s 7 b(t*). 
The increase in H due to the extra energy density will 
shift in the CMB peaks to smaller angular scales (higher I 
values) and increase the peak separation. The extra radi- 
ation also lowers the redshift of matter-radiation equality 
which will increase the suppression of the growth of CDM 
perturbations on subhorizon scales prior to matter domi- 
nation and cause an increase in the temperature variance 
[51]. The increase in the Hubble parameter also reduces 
the contribution to the anisotropy spectrum from the ve- 
locity perturbation and the Silk damping is also less ef- 
fective, further increasing the peak enhancement at large 
I. Lastly, since the extra energy density associated with 
"A" can cluster under the influence of gravity, increasing 
AN U enhances the early Integrated Sachs- Wolfe (ISW) 
effect, leading to a significant increase in the height of 
the first peak relative to the others. 

"Q" 

In general the perturbations in the Quintessence en- 
ergy density, pressure, and peculiar velocity, after Fourier 
transforming to k-space are 

SpQ = ^6^ + V4,6(f>, (39) 

SPq = ^ty-fyty, (40) 
v Q = k6<f>/j> (41) 

and the Klein-Gordon equation for the evolution of the 
perturbations is simply 

54> + 2H5<j)+ (k 2 + a 2 V 4><p )5<t>-lj)5 c = 0. (42) 

There is no shear term for the Quintessence field. The 
density perturbations in the field are not adiabatic 
[52, 53, 54]. That is, 5P Q - c 2 Q 6p Q ^ where c 2 Q = 
dPQ/dpQ is the adiabatic sound speed, a fact which, ac- 
cording to Ratra & Peebles [54] is "exceedingly fortu- 
nate" since pressure fluctuations can resist the collapse 
of Quintessence density fluctuations even if the equation 
of state is zero. In principle the CMB anisotropy spec- 
trum emerging from a quintessence model is very dif- 
ferent from the case of "A" , but certain effects may be 



missing or small depending on the exact behavior of the 
field. 

The changes wrought upon the CMB anisotropy spec- 
trum from replacement of the vacuum energy (A) with 
a Quintessence dark energy that follows the exponen- 
tial potential were investigated by Ratra & Peebles [54], 
Ferreira & Joyce [6] and Skordis & Albrccht [52]. Deep 
in the radiation dominated epoch the AS potential may 
be approximated by a pure exponential and for this po- 
tential, on supcrhorizon scales and with adiabatic initial 
conditions, there is no change to the evolution of per- 
turbations, i.e., 5c ot k 2 t 2 [6], which is wholly expected 
because during this period the Quintessence equation of 
state is exactly the same as that of the radiation. The 
adiabatic initial density fluctuation in Q is simply [6, 52] 

Sq = ^5 c (43) 

which translates to a perturbation in the field [52] 

Scj) = 4<5c/5A 6<j> = 4 8 c /bX. (44) 

Though we have used equations (43) and (44) in all our 
calculations the evolution of the Quintessence perturba- 
tions is largely independent of the exact initial conditions 
of the field when the equation of state is a constant [55] . 

Life becomes more interesting on subhorizon scales, 
kr 3> 1. As in Ferreira & Joyce [6], we take the potential 
derivative and gravitational feedback terms in equation 
(42) to be negligible relative to the k 2 Sq. piece so that, 
after using 3(wf + 1) — 2 = 2/Ht, we obtain 

r2 ^ + 3 (% 4 + T l)-2 ^ + fc2T2 ^"°- (45) 

The solutions of this equation are linear combinations 
of J p (/ct)/t p , N p (kT)/T p with p = 1/2 during radiation 
domination and p = 3/2 during matter domination [6]. 
The subhorizon perturbations in the field oscillate with 
decaying amplitudes and so Quintessence does not clus- 
ter. The density contrast of the CDM continues to grow 
as In t during radiation domination, but during matter 
domination the evolution changes to 5c oc r 2+e where 

2e = 5yjl - 24 n^ D) /25 - 5 < 0. (46) 

The lack of clustering in the Quintessence energy density, 
along with its contribution to H inhibit growth in the 
CDM. 

These results apply when the Albrccht- Skordis poten- 
tial may be approximated by a pure exponential but 
when the polynomial prefactor becomes important, these 
results break down. However we can anticipate some of 
the effects of the field from the simple fact that the effect 
of a significant Quintessence energy density is to drive a 
swifter expansion. In this regard the "Q" model we have 
adopted bears some similarity to "A" in that the effect of 
the extra energy density in the early universe leads to a 
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smaller sound horizon s 7 b(t*) at recombination driving 
the anisotropy spectrum peaks to smaller scales and in- 
creasing their separation. At the same time the increase 
in H also leads to a decrease in the Silk damping, thus 
increasing the temperature variance at large t relative to 
the first peak. However, unlike "A" , Quintessence may 
also significantly reduce the conformal angular diameter 
distance to last scattering thus partially mitigating the 
shift in the location of the first CMB peak [50] . For the 
AS potential, this is not expected to be a large effect be- 
cause the field becomes trapped in its minimum at z <~ z\ 
which, for ^q ?D ' ) 0.1, is above z\ > 16. Another major 
difference between "A" and "Q" is the absence of cluster- 
ing on subhorizon scales, hence there is no enhanced early 
ISW effect for Quintessence. But there is an ISW effect 
for this potential [52] because of the dramatic decrease 
in pq as the minimum of the Quintessence potential is 



approached before ~ za ■ This results in an intermediate- 
to-late ISW effect reducing all the peak amplitudes in the 
power spectrum because the COBE normalization fixes 
the amplitude of the power spectrum at large scales. 

"G" 

Here, once again, there are two different approaches 
to perturbations in the field depending upon the prefer- 
ence for unmodified forms of the geometric/gravitational 
quantities or for unmodified stress-energy terms. As be- 
fore, we adopt the latter approach. From this vantage 
point the perturbations in the energy density, pressure, 
and velocity divergence are the same as in equations (39), 
(40), and (41), but the Klein-Gordon equation for the 
perturbations now assumes the form 



2 H 54> + (fc 2 + a 2 V^) 64>-4>6 c 
a 2 R 



2E 



F 2 3jffi 
00 F 2F 



a 

y 

a 2 F 



2FE 



[F+6R + FmR5<I>] 

Sp - 35P + 3(F + F u %) , 



<t> 



6F, 



(47) 



term. Note that for this 



6F A 



and we have dropped the F^, 
potential and coupling 

Z{F^V U + F^V^Stj)- 

11 

= -3£ (1 + a) 6p G -3£(1- a) 6P G . (48) 

The minimally-coupled limit of this model is another fre- 
quently studied quintessence potential [54, 56, 57]. This 
potential is in some respects the opposite of the Albrecht- 
Skordis model where the energy density during the early 
Universe can be considerable and may lead to a swifter 
expansion prior to its entrapment at z\. In contrast, 
for the Inverse Power Law potential the opposite occurs: 
the Quintessence energy density is inconsequential dur- 
ing much of the evolution of the Universe, only becoming 
important as the present epoch is approached. 

In the minimal-coupling limit the supcrhorizon pertur- 
bations during the radiation dominated epoch may be 
derived in the usual way, leading to 



Sg 



4a 



3 (5a + 26) 



5a + 26 



5c 



8 + 2a 
5a + 26 



5c 



(49) 



(50) 



with no change in 5c- Note that when a^cowe re- 
gain the results of the exponential potential (see equa- 
tion (43)) and, when a = 0, the vacuum result 5\ = 0. 



Again, on supcrhorizon scales there is no change in the 
evolution of 6c during either matter or radiation domina- 
tion and similarly, on subhorizon scales, the Quintessence 
perturbations decay [6, 54] but unlike U Q" there is no 
suppression in the growth of matter perturbations. The 
negligible Quintessence energy density also means that 
there is no change in the size of the sound horizon at 
recombination. There is, however, a significant decrease 
in the angular diameter distance because of the contri- 
bution from the Quintessence energy density at low red- 
shifts that causes a shift in the peak positions to smaller 
I and decreases their separation. This effect increases 
with a because the equation of state and, thus, the en- 
ergy density pc are correlated with the exponent: smaller 
a correspond to smaller pc at a fixed redshift. The peak 
heights are suppressed relative to the COBE normaliza- 
tion point because the larger Quintessence energy density 
terminates matter domination at an earlier epoch. This 
enhances the late ISW effect, increasing the variance on 
large scales and therefore lowering the initial amplitude 
of the power spectrum. Once again, this effect increases 
with a for exactly the same reason. 

For the non-minimally coupled case, £ ^ 0, the situ- 
ation is much more complex. All the effects discussed 
above for the minimally-coupled limit apply but now we 
must also take into account the change in F. A reduction 
in the strength of gravity, F > 1, now slows the expan- 
sion, leading to a larger sound horizon at recombination. 
If F remained constant there would be a compensating 
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change in the distance to the last scattering surface and 
the peak positions would be unaffected. However F must 
decrease in order to attain F = 1 at the present time, so 
this distance can only be smaller than that required to 
leave the peak positions unaltered. Hence a decrease in 
G shifts the CMB anisotropy spectrum peaks to even 
smaller I. The changing F also reduces the ratio of the 
amplitude of the first peak to the others by partially can- 
celling the early ISW effect. Simultaneously, the decrease 
in H shifts the peaks to larger scales due to the fact that 
thermal contact between the baryon-photon fluid can be 
maintained down to a lower the redshift. The slower 
expansion also increases the duration of recombination 
leading to more effective photon diffusion, decreasing the 
temperature variance on small scales. 



VI. CONSTRAINTS FROM THE CMB 

Unlike their similarities for BBN, the three cases we 
are considering here, extra relativistic energy density 
("X"), a minimally coupled scalar field ("Q"), and a 
non-minimally coupled scalar field ("G"), all influence 
the formation of the CMB anisotropy spectrum differ- 
ently, albeit with some features in common. To explore 
the CMB constraints we use modified versions of CMB- 
FAST [58] to construct the anisotropy spectra and vary 
AN„ and the baryon density parameter 7710- We com- 
pare the models to the data from the BOOMERANG 
[59], MAXIMA [60], and DASI [61] detectors, employing 
RADPACK [62] to determine the goodness of fit and we 
assign the confidence level based on Ax 2 - RADPACK 
also allows us to adjust the calibration of each data set, 
at the cost of a x 2 penalty. The current data covers the 
first three peaks in the Cg with reasonable accuracy, so 
while subtle effects may be missed, this approach is suf- 
ficient to extract the gross features of each model. Since 
here we are concentrating on constraints in the AN V - 
77 plane, we have limited the priors, adopting those of 
the "concordance ACDM" model: 51 to t = 1, = 0.3, 
H = 72 kms^Mpc" 1 , along with no "tilt" (n = 1). As 
in related, earlier work (see, e.g., [63, 64, 65, 66]) here we 
explore the constraints on AA„ derived from the CMB 
anisotropy. 

Because the CMB constraints in the non-minimally 
coupled case "G" differ from those for "A" and "Q" we 
discuss this case, G, separately. We note here that Chen 
and Kamionkowski [71] have discussed how constraints 
from the CMB may be used to provide a test of Brans- 
Dicke cosmology, which bears some relation to our case 
"G". 

First, we consider "A". In general, there is an inti- 
mate interplay between 7710 and AN* as shown by the 
solid curves in Figure 4 where our results reveal that an 
increase in AN V can, to some extent, be compensated by 
a reduction in n. This degeneracy is the result of similar 
effects upon the relative heights of the first and second 
peak and is only broken by the presence of the third peak 



in the spectrum. Increasing AN* leads to a larger early 
ISW effect boosting the variance at the first peak rela- 
tive to the rest, while increasing 7710 changes the relative 
heights of the odd and even peaks. At the same time 
increasing AA^- reduces the sound horizon by reducing 
the conformal age of the Universe at last scattering, but 
the presence of pB in equation (38) reveals that a re- 
duction in 7710 can compensate this change by increasing 
the baryon-photon sound speed which overwhelms the si- 
multaneous increase in redshift at which recombination 
occurs. 




FIG. 4: CMB likelihood contours, 68%, 95%, and 99%, for 
"A" (solid) and "Q" (dashed). The best fit points are indi- 
cated by the cross for "X" , and by the square for "Q" . 



The dotted curves in Figure 4 show the contours for 
"Q" . In this case they lack the strong anticorrelation 
between AN„ and 77 that "A" exhibited. Despite the 
differences between "A" and "Q" , the bounds from the 
CMB are quite similar in both cases, with best fit values 
of r/io w 6.2 - 6.5 (n B h 2 w 0.023 - 0.024) and non-zero 
AA„ w 1.2 - 1.6, while the 99% upper bound to AA„ in 
both cases is AA„ < 3.2 (Q { q D) < 0.34, A > 3.4). 

For comparison with the CMB fluctuation data, the 
case of non-minimally coupled fields has several unique 
features which distinguish it from the other two cases 
considered above. The confidence contours shown in Fig- 
ure 5 compare the results for two choices of a. A compar- 
ison of Figures 4 & 5 reveals that the anticorrelation be- 
tween 7710 and AN^ is very strong for the non-minimally 
coupled case G, with the centroid of the contours depen- 
dent on a as is the elongation of the contours. We also 
note that the % 2 at the best fit point decreases with a 
i.e., a — 2 is a better fit than a = 4 and a = 1 is a bet- 
ter fit than a = 2. The large shift in and the sensitivity 
of 7710 to a can be understood as follows. Inverse Power 
Law potential models gcnerically reduce the heights of 
all the peaks because of the late ISW effect and this sup- 
pression decreases with a. The peak suppression favors 
shifts in the experimental calibrations, so by increasing 
the baryon to photon ratio there is a deflection of some 
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of the \ 2 penalty by improving the fit to the first and 
third peaks at the expense of making the fit to the sec- 
ond peak worse. Additionally raising 7710 suppresses the 
Silk damping boosting the temperature variance at large 
t relative to that at the first peak. Lastly as we noted 
before for "A", raising 7710 will also reduce the sound 
horizon at last scattering which will compensate for the 
reduction of the angular diameter distance in these mod- 
els. The shift in the centroid of the contours is due to the 
reduction of the strength of the late ISW effect as a de- 
creases since the increase in 7710 needed to raise the peak 
amplitudes becomes smaller. The anticorrelation of 7710 
and AN? is largely the result of the change in the early 
ISW effect. As AN? increases the ratio of the first to 
second peak amplitudes also increases which, like AN? , 
can be compensated by reducing 7710- 
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FIG. 5: CMB likelihood contours, 68%, 95%, and 99%, for 
"G" and two values of a, the IPL exponent. The solid con- 
tours are for a = 4 and the dashed contours are for a — 1. 
The best fit points are indicated by the diamond for a = 4 
and the cross for a = 1. 



Clearly, in this case the best fit values and the bounds 
for 7710 and AN? are functions of a. For a — A we have 
7710 = 10.9, AN? = -0.84 with a 99% upper bound of 
AN? < 1.8, for a = 2 we find 7710 = 10.0, AN? = -0.90 
with a 99% upper bound of AN? < 0.2 while for a = 1 
we find r/10 = 8.8, AN? = -0.63 with a 99% upper 
bound of AN? < 0.3. However, as we mentioned earlier, 
the \ 2 a t the minima for a = 4 is large (130 for 39 
degrees of freedom) and we can rule out this model as 
incompatible with the data an extremely high confidence 
level (formally, at <~ 10~ n ). The case of a = 2 is only 
marginally CMB-compatiblc, at the 99% confidence level, 
while the a = 1 model is perfectly acceptable. These 
results are in agreement with previous studies on limits 
to a: Malquarti and Liddlc [67] find a < 2 at 95% in the 
minimal limit while Bean & Melchiorri [68], Hannestad 
& Mortsell [69] and Melchiorri et al. [70] , have found that 
the effective equation of state for Quintessence must be 
close to —1 which imposes a similar constraint on a [19]. 



As in our discussion of BBN, we must discuss some 
caveats to our CMB results. When AN? > 0, the en- 
hanced relativistic energy density keeps the early uni- 
verse radiation-dominated to a lower rcdshift. To some 
extent this can be compensated by an increase in the 
matter density Qm- However, since in our CMB fits we 
have fixed at 0.3, the ranges of AN* shown in Fig- 
ure 4 are overly restrictive (compare, e.g., with [66]). 
However, this bias will be ameliorated in the combined 
BBN - CMB likelihood distributions since, as may be 
seen by comparing Figures 3 and 4, the dominant con- 
straint on AN? is provided by BBN. We must again 
point out that the solar system limits to G' will further 
restrict the allowed range of AN? and indeed the best 
fit points shown in Figure 5 are all outside of the allowed 
range corresponding to the bounds in G in eq. (19). For 
non-minimally coupled scalar fields the upper and lower 
bounds to AN? are functions of the Inverse Power Law 
exponent a in the adopted scalar field potential. 



VII. COMBINING BBN AND CMB 

For a self-consistent cosmology any modifications to 
SBBN at rcdshift <~ 10 10 must be consistent with those 
deviations from the standard-model predictions for the 
later evolution of the universe as probed by the CMB 
(10 4 > z > 0). While BBN and the CMB favor slightly 
different regions in the AN„ - 77 plane, there is, in- 
deed, overlap for models with new particles ("A") and 
for those with minimally-coupled scalar fields ("Q"). In 
Figure 6 are shown the joint BBN-CMB likelihood con- 
tours for these two cases. Although the best fits do 
occur for small, nonzero values of AN V (AN U w 0.1, 
7710 w 6.6 <^=> fl B h 2 w 0.024), the deviation from AN V = 
is not statistically significant and the "standard" model 
{AN V = 0) is entirely consistent with the BBN and CMB 
data. For both "A" and "Q", the minimum x 2 in our 
combined BBN - CMB fit is 47 for 39 degrees of freedom; 
for these cases there is consistency between the models 
and the universe at z w 10 10 and at z w 0. 

In top panel of Figure 7 are shown the joint BBN- 
CMB likelihood contours for the "G" case, but only for 
a = 1: we do not show likelihood contours for the other 
two a values since these models (a = 2, 4) are ruled 
out at better than 99% confidence. For a = 1 the y? 
at the minimum for the joint BBN-CMB fits is slightly 
larger than for "A" and "Q", at 57 for 39 degrees of 
freedom; this non-minimally coupled model is compatible 
with the data at 97% confidence level. The minimum for 
this model occurs at AN? w -0.51, 7710 ~ 8.2 (Q B h 2 w 
0.030), with a significant autocorrelation between the two 
quantities. The deviation from AA„ = is now rather 
large: the minimum along AA^ — is located on the 
90% confidence contour. 

Finally in the bottom panel of Figure 7 we impose 
the limits to G' we have mentioned so frequently. The 
minimum from the top panel is outside the allowed region 
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FIG. 6: Joint BBN and CMB likelihood contours, 68%, 95%, 
and 99%, for "X" (solid) and "Q" (dashed). The best fit 
points are indicated by the cross for "X" and by the square 
for "Q". 



AN? = ±0.15 and now resides on the boundary at AN? 
w -0.15, 7710 « 7.6 (fi B /i 2 « 0.028). The x 2 at the 
minimum is such that the goodness-of-fit of the model is 
only compatible with the data at the 98% level. 



VIII. SUMMARY 

BBN and the CMB provide complementary constraints 
on models for the dark matter/energy at completely dis- 
tinct epochs in the evolution of the universe. Consistency 
with both sets of constraints can help in distinguishing 
among different models for new physics beyond the stan- 
dard models of cosmology and particle physics. As dis- 
cussed here, models with extra particles which are rela- 
tivistic during the early evolution of the universe ("A") 
and those with minimally-coupled scalar fields ("Q") 
yield similar predictions for BBN (7710 w 5.0 and AA„ w 
0; see §111, especially Figs. 1 & 3) and for the CMB 
(7710 « 6.2-6.5 and AA„ w 1.2-1.6; see § VI and Fig. 4). 
As a result, they also provide similar good fits to the 
joint BBN-CMB constraints (see §VII and Fig. 6) with 
the baryon density being largely determined by the CMB 
constraints and AN U by BBN. While the best fits to the 
data for these two options are in fact compatible with 
AA„ = (at flBh 2 w 0.024 in both cases), some dark en- 
ergy is permitted: at 99% confidence level AN?< 1.07 
while AN?< 0.85 («[f D) < 0.12, A> 5.7). In contrast, 
the new physics associated with non-minimally coupled 
scalar fields ("G") induces some unique behavior, leading 
to quite different BBN and CMB constraints which, in 
general, are quite difficult to satisfy simultaneously Be- 
cause such models generally modify the strength of grav- 
ity during the earlier evolution of the universe, A^T < 
is not only allowed, it is favored by both the BBN and 
CMB constraints (see §111, especially Fig. 2 and §VI, es- 
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FIG. 7: Joint BBN and CMB likelihood contours, 68%, 95%, 
and 99%, for "G", a = 1. The top panel is with no G con- 
straints eq. (19), the bottom is with those constraints. The 
best fit points are indicated by the crosses. 



pecially Fig. 5) . However, for "G" the allowed regions in 
the AN V - 77 plane are sensitive to the form of the scalar 
field potential (see Fig. 5). For inverse power law poten- 
tials, only those models with a< 1 have CMB-idcntified 
regions which have significant overlap with the regions 
compatible with the BBN constraints (compare Figs. 2 & 
5). Furthermore, since there are solar system constraints 
on the possible variation of G, there are non-BBN and 
non-CMB constraints on G'/G (see eq. (19)) which pro- 
vide independent constraints on AN?. When these are 
combined (for a = 1) with the BBN and CMB constraints 
the best fit value of AN? is at -0.15 and Q B h 2 w 0.028 
(see Fig. 7). However, we note that there is less than 
a 2% probability that this model is, in fact, compatible 
with the current BBN and CMB data. 

The two models for Quintessence we have investigated 
in this paper are by no means the only plausible exam- 
ples of dynamic dark energy. What distinguishes the 
two cases we have considered is their importance dur- 
ing BBN and their role in the generation of the primary 
anisotropies in the CMB. In fact, in this regard, these 
two models are unique since it is much more common for 
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Quintessence models to become important only at low 
rcdshifts. In such cases BBN provides few, if any, con- 
straints and the primary anisotropies of the CMB are 
unaffected. The influence of the dark energy is via the 
change in redshift of matter-dark energy equality and 
thus the corresponding COBE normalization of the mat- 
ter power spectrum. For the CMB this translates into 
a global stretching of the peaks in the spectrum and a 
reduction in their amplitudes but for Large (and Small) 
Scale Structure (LSS), Type la supernovae observations, 
and weak lensing surveys the change in the extent of 
matter perturbation growth is more important and their 
observation will provide valuable additional and compli- 
mentary constraints upon the models [72, 73, 74]. The 
minimally coupled IPL potential falls into this category. 
The CMB already constrains a < 2 [67] and including 
LSS and Type la supernova data yields wq ^ — 0.7 at 
the present time which translates into a < 1.5 for the 
concordance model values of Ho and Om- 

These same cosmological tests also furnish constraints 
upon the two models we have examined: "Q" and"G". 
For the pure exponential potential these signatures were 
discussed by Ferreira & Joyce [6] and the Albrecht Sko- 
rdis modification (our case "Q") by Skordis & Albrect 



[52]. In contrast with typical Quintessence models the 
AS potential modifies the matter power spectrum mainly 
by inhibiting the growth of density perturbations rather 
than changing the COBE normalization. The essentially 
stationary field below the redshift of 15 (correspond- 
ing to A> 5.7) leads to little change in the magnitudc- 
redshift relation that could be probed by Type la super- 
novae. For the non-minimally coupled Inverse Power Law 
model the effects are supplementary to those in the min- 
imal limit. A non-zero £ will introduce additional effects 
through the change in gravitational strength F that in 
turn will alter the growth of density perturbations by di- 
rectly influencing their growth in addition to the change 
in the normalization. 
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